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ABSTRACT 

We discuss the possibility of detecting the presence of primordial black holes (PBHs), such as those that 
might account for galactic dark matter, using modification of pulsar timing residuals when PBHs pass within 
^1000 AU and impart impulse accelerations to the Earth. With this technique, PBHs with masses around 
MpBH ~ lO^^g 0.1 lunar mass) can be detected. Currently, the constraints on the abundance of such dark 
matter candidates are weak. A 30 year-long monitoring campaign with the proposed Square Kilometer Array 
(SKA) can rule out a PBH fraction more than ^ 1 / 10 in the solar neighborhood in the form of dark matter with 
massMpBH ^ IQ^^g. 

Subject headings: pulsars: general — black hole physics — dark matter 



1. INTRODUCTION 

Despite a long and intense observational quest, the physi- 
cal identity of dark matter that accounts for most of the mass 
in our galaxy is still unknown. Existing candidates for dark 
matter can be found over many orders of magnitude in mass 
from particle-physics scales to physically large astronomical 
scales (see e.g. Bertone et al. 2005, Carr 2005 for reviews 
on candidates). Among these, primordial black holes (PBHs) 
provide an important astrophysical candidate for dark matter. 
Their abundance is now constrained relatively well for masses 
above Mpbh ^ 10~^g using the galactic microlensing optical 
depth (Alcock et al. 1998) and for Mpbh < lO'^g using the 
Galactic gamma-ray background intensity compared to the 
expected intensity of Hawking radiation associated with the 
black hole evaporation (Carr & Sakellariadou 1999). 

In the mass range lO'^g <Mpbh S 10^"g the femto-lensing 
technique can be used to monitor the presence of PBHs, de- 
pending on structure of gamma-ray-bursts (Gould 1992). Un- 
fortunately, in the mass range lO^^g < Mpbh ^ lO^^g (or 
around 10"^ to 10 lunar masses), there is currently no pre- 
ferred method to constrain PBH abundance (Carr & Sakellar- 
iadou 1999). At these large masses, the flux of PBHs around 
the solar neighborhood is much smaller than the presumed 
flux of particle physics dark matter candidates. Therefore, 
direct detection of PBHs requires detectors with collecting ar- 
eas or cross-sections that span planetary distances of order an 
astronomical unit or larger In such a scenario, the presence 
of PBHs is inferred through gravitational interaction between 
passing PBHs and "test masses" that can sense impulse accel- 
erations. With the typical distances expected for PBHs, these 
impulse accelerations have very small amplitudes and can be 
detected only with high-precision monitoring of test mass lo- 
cations. This approach, however, does not involve compli- 
cated physical processes nor detailed astronomical models 
and enables us to study the relevant parameter space with sim- 
ple physics. 

Accurate monitoring of the location of a set of test masses 
either from ground or space is now underway (or planned) for 
gravitational wave studies with small strain amplitudes (e.g. 
Cutler & Thorne 2002, see also Nakamura et al. 1997, loka 
et al. 1998, Inoue & Tanaka 2004 for detecting gravitational 
waves from PBH binaries). It is possible to detect PBH signa- 
tures with gravitational wave interferometers (Seto & Cooray 
2004, Adams & Bloom2004). Given the time scale of the 



passage and the fact that ground-based detectors are severely 
affected by seismic noise below 10 Hz, these observations re- 
quire space-based gravitational wave detectors, such as the 
Laser Interferometer Space Antenna (LISA) (Bender et al. 
1998). Assuming that the designed acceleration noise level 
around 10~^Hz can be extrapolated down to ^ 10"^ Hz, the 
expected detection rate for LISA is 0.01 events per year for a 
PBH mass of lO'^g. While this detection rate with LISA is 
low, second-generation space detectors might provide inter- 
esting limits (Seto & Cooray 2004). 

In this paper, we discuss an additional technique to detect 
PBHs focusing on pulsar timing residuals. With pulsar timing 
observations, we can measure the pulse-like acceleration by a 
PBH over a time scale up to ^10 yr (or equivalently down to 
3nHz in frequency space). Furthermore, with a Pulsar Tim- 
ing Array (PTA) that monitors a large number of pulsars to 
search for correlated signals buried in the noises of individ- 
ual pulsars ( Hellings & Downs 1983, Foster & Backer 1990, 
Cordes et al. 2004, Kramer et al. 2004, Jenet et al. 2005), we 
can improve the overall sensitivity to detect PBH signatures. 
In this respect, we discuss the feasibility of using proposed 
projects such as the Parkes Pulsar Timing Array (PPTA') and 
the Square Kilometer Array (SKA") to study what constraints 
they place on the fraction of PBHs that make up the dark mat- 
ter While we limit our discussion in this paper to dark matter 
candidates in the form of PBHs, our calculations equally ap- 
ply for other dark matter candidates that are sufficiently com- 
pact and have masses that are similar to the values we study 
here. 



2. EXPECTED DETECTION RATE 

First, we study the passage of a PBH with mass Mpbh, ve- 
locity V, and closest approach distance D from the Sun. A 
schematic diagram of the interaction is given in Figure 1. 
Hereafter, the velocity V is fixed at V = 350km sec"' which 
is the typical velocity for halo dark matter relative to the solar 
system (Carr & Sakellariadou 1999). Besides the combination 
of parameters (Mpbh,D), we can also characterize the event 
with the time scale and the amplitude of the gravitational per- 
turbation induced by the PBH. 

' See http://www.atnf.csiro.aii/resea rch/pulsar/psrtimel 
^ See http://www.skatelescope.org 
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The time scale of the gravitational interaction is 

T = D/V= 13.6yr (D/IOOOAU) (y/350km sec"')"' . (1) 

As we shall see later, the detectable distance range D from 
the Sun for the upcoming pulsar timing arrays is around 1000 
AU. For a circular orbit with radius r around the Sun, the or- 
bital period is 170(r/30AU)-'/^yr and the orbital velocity is 
5(r/30AU)-i/2km sec ' . Note that 30 AU is the distance be- 
tween Neptune and the Sun. Therefore, both the gravitational 
focusing of a PBH by the Sun is negligible and the time scale 
of gravitational perturbation by a bound object at distance 
> lOOAU is much longer than a realistic expected observa- 
tional period which can be as long as several decades (see 
also Zakamska & Tremaine 2005). 

We summarize our parameters and results in Figure 2. 
In this figure the time scale for interaction T is shown 
as the right-hand vertical axis taking the velocity V to be 
350km sec"'. We also show two observational periods with 
Tobs = 1 and 10 yr using horizontal dashed lines. As the rel- 
evant distance D from the Sun is much larger than those be- 
tween planets and the Sun, the acceleration by the PBH should 
be regarded as that for the solar system bary center (SSB). The 
tidal acceleration of the Earth relative to the SSB is suppressed 
by a factor (lAU/D) and is further reduced when integrating 
many orbital cycles. Therefore, we focus on detecting the 
bulk acceleration of the SSB by a passing PBH using pulsar 
timing residuals (see e.g. Damour & Taylor 1991 for a discus- 
sion on the constant SSB acceleration). 

With the time scale established, we now discuss the am- 
plitude of the gravitational perturbation induced by a pass- 
ing PBH which, in turn, leads to a signal in Time-of- Arrival 
(TOA) residuals of pulsars. The amplitude of the pulse-like 
acceleration of a passing PBH with mass Mpbh, distance D, 
and velocity V is GM^yhi/D^, with the duration given by 
=D/y. Thus, the TOA residual tr induced by a PBH is 

5x IGMpBH 2 GMmy^ 
c c D-^ cV^ 

where the quantity 5x is the modulation of the SSB's position 
by the passage of a PBH. For a given velocity V, we have a 
one-to-one correspondence between the mass of a PBH Mpbh 
and the pulsar timing residual f,.. Using appropriate numbers, 
the timing residual is 

In Figure 2 the residual tr is indicated by the top horizon- 
tal axis on the mass-distance plane, and again taking V = 
350km sec"'. From the observed time scale Tobs and the tim- 
ing residual t^, we can inversely estimate the distance D and 
the mass Mpbh, assuming a magnitude for the velocity V . 
Note that, in a general situation, there are two observable 
quantities with three free parameters. It is difficult to deter- 
mine all three quantities D, Mpbh and V simultaneously by 
observing a pulse-like acceleration, unless one of these pa- 
rameters can be established independently. Here we fix the 
velocity, but with any detection, this degeneracy can impact a 
simple interpretation of an event. 

Next we estimate the flux of nearby PBHs. The den- 
sity ppbh of PBH around the Sun must be smaller than 
the estimated local density of dark matter pdm = 0.011 ± 
0.005[Mopc-3] (OlUng & Memfield 2001). We denote the 
fraction of dark matter in the form PBHs as /pbh = Ppbh/pdm- 



For a given mass Mpbh, the event rate R of PBH passing near 
the Sun within a distance D is R = ttppbhD^V /Mpbh, or 

/? = 0.068yr-' ( ^) ( EE^) 

^ \l.O J VO.OIIM0PC-3 J 

( y \ /MpbhVY d y 

^ ^350kmsec-' / \ \(fl^g) VlOOOAUy 
In Figure 2 we show the distance D for the given event rates 
R = Q.\ yr"' and lyr"' under the assumption that /pbh = 1 
(Ppbh = Pdm)- Here, note that the distance depends on the 
mass as D oc Mpg"^. 

If our observational time period for the monitoring of pul- 
sars Tobs is less than 10 yr, we need to search for events 
below the upper short-dashed line in Figure 2. At the same 
time, for an event to be present given the expected event rate, 
it should be above the lower long-dashed line. The result- 
ing mass range is Mpbh ^ 4 x lO^^g corresponding to a TOA 
residual of f, < 7ns. This is much smaller than the expected 
noise level ^ 100ns for relatively stable pulsars with a 10 yr 
observational period {e.g. Jenet et al. 2005). Therefore, it 
is statistically difficult to detect a passing PBH using a single 
pulsar. The detection might become possible with a pulsar 
timing array (PTA) where one can analyze a large number of 
pulsars to search for a correlated timing residuals below the 
noise level of individual pulsar. 

We now discuss the possibility of using PTA projects for 
a detection of PBHs. Pulsar timing observations have been 
long used to search for the low-frequency gravitational wave 
background. The basic theoretical framework is well devel- 
oped (see Detweiler 1979, Hellings & Downs 1983, Bland- 
ford et al. 1984, Foster & Backer 1990, Damour & Taylor 
1991, Kaspi et al. 1994, Thorsett & Dewey 1996, McHugh et 
al. 1996). A pulse-like gravitational wave with a time scale 
T and a non-dimensional amplitude he modulates radio sig- 
nals from a pulsar by /S.v/v ^ he {v: frequency of the radio 
pulse). The same nondimensional observable Sc = Av/v gen- 
erated by the PBH acceleration is given as 

5v 5v 1 GMpBH ^ GMpBH 

sc= --2—T -— , (5) 

u c c cDV 

where the quantity 5v is the velocity modulation of the SSB 

by the PBH. This expression is evaluated as Sc = 4.3 x 

10-'^ (^) (ti^)"' (3511^)"' • In Figi^'-e 2, ampli- 
tude Sc is shown with solid lines. Here, we discuss the modu- 
lation of pulsars' signals over a given time scale T . For refer- 
ence, note that TOA residual f,- and the nondimensional ampli- 
tude Sc are almost equivalent in Fourier space where we have 
t,-- f = Sc with / ~ as seen with eqs.® and (|5]). 

For each pulsar, the TOA modulation due to the gravita- 
tional wave background is given by two terms, the pulsar term 
and the Earth term. The magnitude of these two contributions 
are comparable. The pulsar term is determined by the gravi- 
tational wave background at the pulsar, while the Earth term 
correlates the TOA residuals for all pulsars (Detweiler 1979, 
see also Jenet et al. 2003). The pulsar term can be regarded as 
noise when searching for the correlated signal by the gravita- 
tional wave background with a PTA. The TOA residual by a 
passing PBH around the Earth is completely correlated among 
pulsars. This is advantageous when searching for a PBH sig- 
nature, especially in the case when the signal is strong enough 
to be detected with a small number of pulsars. But this does 
not hold when noise of each pulsar dominates the correlated 
signal and we need a PTA to detect it. 
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Fig. 1 . — Schematic picture for PBH searcli with pulsar timing observation. 
A PBH with mass Mpbh and velocity V passes near the Sun with closest 
distance D 2> lAU. We fix the velocity V at its typical value 350km sec"' . 

Here we calculate the expected sensitivity tr as a function 
of the signal duration T , or equivalently, we examine the re- 
gion in mass-distance plane (Figure 2) that can be probed with 
upcoming projects. To estimate the measurement sensitivity 
tr{T) as a function of the time scale T , we first summarize re- 
sults for the sensitivity h^T) given for an observation of the 
gravitational wave background with a PTA. This will give us 
an estimation of the measurement sensitivity for a PBH event 
in terms of Sc(T). Then, we can derive the sensitivity in terms 
of f,(r) using the simple correspondence between Sc and t,- 
mentioned before. In the next section, we will discuss how 
we can separate PBH signatures from signatures related to a 
gravitational wave background. 

Because a PBH signal has a single pulse-like profile with 
a finite duration T, one might expect that the detection sig- 
nificance for the event is not improved with an observational 
period Tobs much longer than the signal duration T . In re- 
ality, parameters of pulsars such as their directions and dis- 
tances can be estimated better with a long observational pe- 
riod and, consequently, timing residuals are less affected by 
errors in these parameters (Blandford et al. 1984). While this 
dependence on the observational period Tobs is important for 
Tobs ^ 3yr, we neglect it with the assumption of a sufficient 
observational period Tobs > 5yr. 

For studies involving the gravitational wave background. 
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to describe the 



we often use the spectrum flcw(f) — p ji^f 
energy density of waves per logarithmic frequency interval 
around a frequency /, when normalized by the critical den- 
sity pc of the universe (see e.g. Allen 1997, Maggiore 2000). 
The characteristic amplitude hdT) is related to ilcwif) at 
T =/-! as flcwif) = Ts^fhciTf, (Hq = 100/zkm sec-'Mpc-i: 
the Hubble parameter) and we can write 



hc(T) = 4.0 X 10" 



ncw(f=T-')h^V^^f T 



10- 
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VlOyr 



(6) 



For a given observational period Tobs, it is difficult to ex- 
tract information on gravitational waves with periods longer 
than observational duration. This is because we cannot sep- 
arate TOA residuals due to gravitational waves from secular 
effects such as the long-term spin evolution of pulsars (Kaspi 
et al. 1994, see also Damour & Taylor 1991). With the pulsar 
timing analysis, the sensitivity of the spectrum ilcwif) gen- 
erally becomes best at the frequency / ~ 1 /T^hs or time scale 
T ~ Tobs- In terms of the characteristic amplitude hdT), the 
expected sensitivity atT = Tobs can be written with the scaling 
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Fig. 2. — Mass-distance relations under various observational conditions 
for PBH search. The short-dashed lines (D oc Mp^^^) show the signal dura- 
tion for the fixed PBH velocity V = 350km sec"' . We have one to one corre- 
spondence between the mass Mpbh (bottom axis) and the residual time t,- (top 
axis). The same is true for the distance D (right axis) and the time scale T (left 
axis). The solid lines (D oc Mpg^^) show the amplitudes Sc of timing residuals 
in the form iwjv. The triangles coirespond to the points (Tohi-^GViC^^^^^) 
for the following three pulsar timing projects, SKA, PPTA and current results. 
The hatched region is detectable event for SKA bounded by the short-dashed 
and dash-dotted lines. The triangle for SKA moves along the dash-dotted 
line (D oc Mpgpj) for different observational period T„\,^. The long-dashed 
lines are distances D for given event rates R and mass Mpbh with putting 
PPBH = PDM = 0.01 lM0pc"'. 

-3/2 

relation hdJohs) oc T^^l (Rajagopal & Romani 1995, Jaffe 
& Backer 2003), with an overall normalization that depends 
on the parameters of the individual PTA. Similarly, we have 
if (T) oc r~^/^ or f, (r) cx r~'/^ for the PBH acceleration with 
T < Tobs- From eqs.([T]i and dU, the above scaling relation 
becomes D oc Mpgjj in the mass-distance plane. 

For the present analysis the sensitivities of various pul- 
sar timing projects can be characterized by the combination 
(Tobs, ilGw{T~i^f)h^) associated with a gravitational wave back- 
ground search. Here we use values {Tohs,^cw(T~f^Jh^) = 
(Syr, 9.1 X 10"'^) for the two parameters for the PPTA fol- 
lowing Jenet et al. (2006) (ilcwh^- 95% detection rate), and 
(lOyr, 10"'^^) for the SKA-PTA from Damour and Vilenkin 
(2005) (^Icwh^- 6-cr detection significance). As an added 
reference, we also take (10yr,6.0 x 10"^) for the current ob- 
servational sensitivity from Jenet et al. (2006). In Figure 
2, we show the corresponding points (Tobs, civ (T^obil'^^) 
filled triangles on the mass-distance plane using eq.((6l). For 
the SKA-PTA project, the detectable event rate is also shown 
as the hatched region. This region is bounded by the two 
lines related to the requirement T < Tabs and another related 
to D ocMpBjj. 

The triangle for the SKA-PTA project with Tobs = 10 yr is 
nearly on the long-dashed line R = O.lyr"'. This means that 
with /pBH = 1 and Mpbh ~ 3 x lO^'^g, we will be able to de- 
tect one PBH event in a 10 yr observation with a duty cy- 
cle C = 1. With an observational period Tobs ^ 30 yr, the 
triangle for the SKA-PTA moves upward along the dash- 
dotted line and is now nearly on the long-dashed line for 
R = lyr"' and mass Mpbh 2 x 10^"^g. This indicates that 
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if the local dark matter is made by PBHs with /pbh = 1 and 
MpBH ~ 2 X 10^'*g, the pulse-like signals for passing PBHs 
with time scales T ^ Tobs ~ 30 yr will highly overlap in data 
streams of pulsar timing. 

The duty cycle C is proportional to D^T/Mpbh T^'/Mpbh 
for events with a time scale T and a mass Mpbh(oc fr)- If a 
point in the detectable zone (the hatched region in Figure 2) 
on the mass-distance plane is expected to have C > 1 as is the 
case for the previous example, there is an appropriate PBH 
mass above which the duty cycle is lower than 1, and events 
will be resolvable for the same time scale T . We can also 
decrease the expected duty cycle by searching shorter time- 
scale events and keeping the target mass Mpbh- In this case 
the new point with C < 1 on the plane might be out of the 
detectable zone after crossing the left boundary (dash-dotted 
line) from above. 

On the other hand, if no PBH signatures are detected with 
the SKA-PTA in 30 yr, we can constrain their abundance frac- 
tion down to /pbh 1/10 around Mpbh '--^ 2 x 10^''g. Of 
course, this is based on an order of magnitude estimate and we 
need detailed analysis to evaluate the exact fraction. While we 
have discussed the SKA-PTA so far, using the PPTA over an 
observational period Tobs of 5 yr, it will be difficult to detect a 
single PBH. With Fobs 20 yr, the PPTA has the potential to 
provide a useful constraint around Mpbh ^ lO^^g. 

3. DISCUSSIONS 

As we have discussed so far, it is quite important to have 
a long observational period Tobs and analyze long-term TOA 
residuals when searching for PBH signatures. For this pur- 
pose, we need stable reference clocks and high-precision 
ephemerides of the solar system. While we do not discuss 
required accuracies of them, those requirements are similar to 
those that are needed for gravitational wave searches (Foster 
& Backer 1990, Kaspi et al. 1994, Hobbs et al. 2006). 

Merging super-massive black hole binaries are considered 
to be the dominant astrophysical source of the gravitational 



wave background around a frequency of ^ O.lyr" (or a time 
scale T ^ 10 yr) (see also Damour & Vilenkin 2005 for the 
gravitational wave background generated by cosmic strings). 
The estimated amplitude using simple assumptions is he ^ 
10-15-5- lO-"" at /= O.lyr-i (Jaffe & Backer 2003, Wyithe & 
Loeb 2003). But depending on model parameters, a smaller 
amplitude he ^ 10" is also predicted (Enoki & Nagashima 
2006). While the estimated amplitude has a large uncertainty 
at present, this gravitational wave background can be an ef- 
fective noise for a PBH search and vice-versa. Therefore, it is 
desirable to observationally distinguish these two signatures 
in TOA residuals. 

One straightforward approach is to investigate spectral in- 
formation of the residuals expected for the PBH accelerations 
and at the same time the gravitational wave background. An- 
other approach is to use the angular pattern of the residu- 
als on the sky. A PBH signature will have a dipole pattern 
(Z = 1) whose direction is determined by the acceleration vec- 
tor. In this respect the dipole component parallel to the eclip- 
tic plane will be more affected by errors in the solar system 
ephemerides than the component perpendicular to the plane. 
On the other hand, the angular pattern due to the gravitational 
wave background does not have a dipole moment but has mul- 
tipole modes starting from the quadrupole (/ = 2) where 75% 
of the angular power is expected to exist (Burke 1975). There- 
fore, by studying the angular pattern of the timing residuals 
we can, in principle, separate the signatures of PBH events 
from gravitational waves. Once PTAs start to monitor suffi- 
cient pulsars, it is clear these topics must be further investi- 
gated. 
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